Genomic prediction for tick resistance in Braford and Hereford cattle.

One of the main animal health problems in tropical and subtropical cattle production is the bovine tick, which causes decreased performance, hide devaluation, increased production costs with acaricide treatments, and transmission of infectious diseases. This study investigated the utility of genomic prediction as a tool to select Braford (BO) and Hereford (HH) cattle resistant to ticks. The accuracy and bias of different methods for direct and blended genomic prediction was assessed using 10,673 tick counts obtained from 3,435 BO and 928 HH cattle belonging to the Delta G Connection breeding program. A subset of 2,803 BO and 652 HH samples were genotyped and 41,045 markers remained after quality control. Log transformed records were adjusted by a pedigree repeatability model to estimate variance components, genetic parameters, and breeding values (EBV) and subsequently used to obtain deregressed EBV. Estimated heritability and repeatability for tick counts were 0.19 ± 0.03 and 0.29 ± 0.01, respectively. Data were split into 5 subsets using k-means and random clustering for cross-validation of genomic predictions. Depending on the method, direct genomic value (DGV) prediction accuracies ranged from 0.35 with Bayes least absolute shrinkage and selection operator (LASSO) to 0.39 with BayesB for k-means clustering and between 0.42 with BayesLASSO and 0.45 with BayesC for random clustering. All genomic methods were superior to pedigree BLUP (PBLUP) accuracies of 0.26 for k-means and 0.29 for random groups, with highest accuracy gains obtained with BayesB (39%) for k-means and BayesC (55%) for random groups. Blending of historical phenotypic and pedigree information by different methods further increased DGV accuracies by values between 0.03 and 0.05 for direct prediction methods. However, highest accuracy was observed with single-step genomic BLUP with values of 0.48 for -means and 0.56, which represent, respectively, 84 and 93% improvement over PBLUP. Observed random clustering cross-validation breed-specific accuracies ranged between 0.29 and 0.36 for HH and between 0.55 and 0.61 for BO, depending on the blending method. These moderately high values for BO demonstrate that genomic predictions could be used as a practical tool to improve genetic resistance to ticks and in the development of resistant lines of this breed. For HH, accuracies are still in the low to moderate side and this breed training population needs to be increased before genomic selection could be reliably applied to improve tick resistance.

A genome-wide scan for selection signatures in Nellore cattle.

Brazilian Nellore cattle (Bos indicus) have been selected for growth traits for over more than four decades. In recent years, reproductive and meat quality traits have become more important because of increasing consumption, exports and consumer demand. The identification of genome regions altered by artificial selection can potentially permit a better understanding of the biology of specific phenotypes that are useful for the development of tools designed to increase selection efficiency. Therefore, the aims of this study were to detect evidence of recent selection signatures in Nellore cattle using extended haplotype homozygosity methodology and BovineHD marker genotypes (>777,000 single nucleotide polymorphisms) as well as to identify corresponding genes underlying these signals. Thirty-one significant regions (P < 0.0001) of possible recent selection signatures were detected, and 19 of these overlapped quantitative trait loci related to reproductive traits, growth, feed efficiency, meat quality, fatty acid profiles and immunity. In addition, 545 genes were identified in regions harboring selection signatures. Within this group, 58 genes were associated with growth, muscle and adipose tissue metabolism, reproductive traits or the immune system. Using relative extended haplotype homozygosity to analyze high-density single nucleotide polymorphism marker data allowed for the identification of regions potentially under artificial selection pressure in the Nellore genome, which might be used to better understand autozygosity and the effects of selection on the Nellore genome.

Genome scan for meat quality traits in Nelore beef cattle.

Meat quality traits are economically important because they affect consumers' acceptance, which, in turn, influences the demand for beef. However, selection to improve meat quality is limited by the small numbers of animals on which meat tenderness can be evaluated due to the cost of performing shear force analysis and the resultant damage to the carcass. Genome wide-association studies for Warner-Bratzler shear force measured at different times of meat aging, backfat thickness, ribeye muscle area, scanning parameters [lightness, redness (a*), and yellowness] to ascertain color characteristics of meat and fat, water-holding capacity, cooking loss (CL), and muscle pH were conducted using genotype data from the Illumina BovineHD BeadChip array to identify quantitative trait loci (QTL) in all phenotyped Nelore cattle. Phenotype count for these animals ranged from 430 to 536 across traits. Meat quality traits in Nelore are controlled by numerous QTL of small effect, except for a small number of large-effect QTL identified for a*fat, CL, and pH. Genomic regions harboring these QTL and the pathways in which the genes from these regions act appear to differ from those identified in taurine cattle for meat quality traits. These results will guide future QTL mapping studies and the development of models for the prediction of genetic merit to implement genomic selection for meat quality in Nelore cattle.

A simple and efficient method for predicting protein-protein interaction sites.

Computational methods for predicting protein-protein interaction sites based on structural data are characterized by an accuracy between 70 and 80%. Some experimental studies indicate that only a fraction of the residues, forming clusters in the center of the interaction site, are energetically important for binding. In addition, the analysis of amino acid composition has shown that residues located in the center of the interaction site can be better discriminated from the residues in other parts of the protein surface. In the present study, we implement a simple method to predict interaction site residues exploiting this fact and show that it achieves a very competitive performance compared to other methods using the same dataset and criteria for performance evaluation (success rate of 82.1%).

The Star STING server: a multiplatform environment for protein structure analysis.

Star STING is the latest version of the STING suite of programs and corresponding database. We report on five important aspects of this package that have acquired some new characteristics, designed to add key advantages to the whole suite: 1) availability for most popular platforms and browsers, 2) introduction of the STING_DB quality assessment, 3) improvement in algorithms for calculation of three STING parameters, 4) introduction of five new STING modules, and 5) expansion of the existing modules. Star STING is freely accessible at: http://sms.cbi.cnptia.embrapa.br/SMS/, http://trantor.bioc.columbia.edu/SMS, http://www.es.embnet.org/SMS/, http://gibk26.bse.kyutech.ac.jp/SMS/ and http://www.ar.embnet.org/SMS.

Predicting enzyme class from protein structure using Bayesian classification

Predicting enzyme class from protein structure parameters is a challenging problem in protein analysis. We developed a method to predict enzyme class that combines the strengths of statistical and data-mining methods. This method has a strong mathematical foundation and is simple to implement, achieving an accuracy of 45%. A comparison with the methods found in the literature designed to predict enzyme class showed that our method outperforms the existing methods.

Building multiple sequence alignments with a flavor of HSSP alignments

Homology-derived secondary structure of proteins (HSSP) is a well-known database of multiple sequence alignments (MSAs) which merges information of protein sequences and their three-dimensional structures. It is available for all proteins whose structure is deposited in the PDB. It is also used by STING and (Java)Protein Dossier to calculate and present relative entropy as a measure of the degree of conservation for each residue of proteins whose structure has been solved and deposited in the PDB. However, if the STING and (Java)Protein Dossier are to provide support for analysis of protein structures modeled in computers or being experimentally solved but not yet deposited in the PDB, then we need a new method for building alignments having a flavor of HSSP alignments (myMSAr). The present study describes a new method and its corresponding databank (SH2QS--database of sequences homologue to the query [structure-having] sequence). Our main interest in making myMSAr was to measure the degree of residue conservation for a given query sequence, regardless of whether it has a corresponding structure deposited in the PDB. In this study, we compare the measurement of residue conservation provided by corresponding alignments produced by HSSP and SH2QS. As a case study, we also present two biologically relevant examples, the first one highlighting the equivalence of analysis of the degree of residue conservation by using HSSP or SH2QS alignments, and the second one presenting the degree of residue conservation for a structure modeled in a computer, which , as a consequence, does not have an alignment reported by HSSP.

PDB-Metrics: a web tool for exploring the PDB contents

PDB-Metrics (http://sms.cbi.cnptia.embrapa.br/SMS/pdb_metrics/index.html) is a component of the Diamond STING suite of programs for the analysis of protein sequence, structure and function. It summarizes the characteristics of the collection of protein structure descriptions deposited in the Protein Data Bank (PDB) and provides a Web interface to search and browse the PDB, using a variety of alternative criteria. PDB-Metrics is a powerful tool for bioinformaticians to examine the data span in the PDB from several perspectives. Although other Web sites offer some similar resources to explore the PDB contents, PDB-Metrics is among those with the most complete set of such facilities, integrated into a single Web site. This program has been developed using SQLite, a C library that provides all the query facilities of a database management system

The Star STING server: a multiplatform environment for protein structure analysis

Star STING is the latest version of the STING suite of programs and corresponding database. We report on five important aspects of this package that have acquired some new characteristics, designed to add key advantages to the whole suite: 1) availability for most popular platforms and browsers, 2) introduction of the STING_DB quality assessment, 3) improvement in algorithms for calculation of three STING parameters, 4) introduction of five new STING modules, and 5) expansion of the existing modules. Star STING is freely accessible at: http://sms.cbi.cnptia.embrapa.br/SMS/, http://trantor.bioc.columbia.edu/SMS, http://www.es.embnet.org/SMS/, http://gibk26.bse.kyutech.ac.jp/SMS/ and http://www.ar.embnet.org/SMS.

STING Contacts: a web-based application for identification and analysis of amino acid contacts within protein structure and across protein interfaces.

UNLABELLED: Amino acid contacts in terms of atomic interactions are essential factors to be considered in the analysis of the structure of a protein and its complexes. Consequently, molecular biologists do require specific tools for the identification and visualization of all such contacts. Graphical contacts (GC) and interface forming residue graphical contacts (IFRgc) presented here, calculate atomic contacts among amino acids based on a table of predefined pairs of the atom types and their distances, and then display them using number of different forms. The inventory of currently listed contact types by GC and IFRgc include hydrogen bonds (in nine different flavors), hydrophobic interactions, charge-charge interactions, aromatic stacking and disulfide bonds. Such extensive catalog of the interactions, representing the forces that govern protein folding, stability and binding, is the key feature of these two applications. GC and IFRgc are part of STING Millennium Suite. AVAILABILITY: http://sms.cbi.cnptia.embrapa.br/SMS, http://trantor.bioc.columbia.edu/SMS, http://mirrors.rcsb.org//SMS, http://www.es.embnet.org/SMS and http://www.ar.embnet.org/SMS (Options: Graphical Contacts and IFR Graphical Contacts).

Defining 3D residue environment in protein structures using SCORPION and FORMIGA.

SUMMARY: Two web-based applications to analyze amino acids three-dimensional (3D) local environment within protein structures-SCORPION and FORMIGA-are presented. SCORPION and FORMIGA produce a graphical presentation for simple statistical data showing the frequency of residue occurrence within a given sphere (defined here as the 3D contacts). The center of that sphere is placed at the Calpha and at the last heavy atom in the side chain of the selected amino acid. Further depth of detail is given in terms of a secondary structure to which the profiled amino acid belongs. Results obtained with those two applications are relevant for estimating the importance of the amino acid 3D local environment for protein folding and stability. Effectively, SCORPION and FORMIGA construct knowledge-based force fields. The difference between SCORPION and FORMIGA is in that the latter operates on protein interfaces, while the former only functions for a single protein chain. Both applications are implemented as stand-alone components of STING Millennium Suite. AVAILABILITY: http://sms.cbi.cnptia.embrapa.br/SMS, http://trantor.bioc.columbia.edu/SMS, http://mirrors.rcsb.org/SMS, http://www.es.embnet.org/SMS and http://www.ar.embnet.org/SMS. [options: Scorpion, Formiga]

ConSSeq: a web-based application for analysis of amino acid conservation based on HSSP database and within context of structure.

SUMMARY: A web-based application to analyze protein amino acids conservation-Consensus Sequence (ConSSeq) is presented. ConSSeq graphically represents information about amino acid conservation based on sequence alignments reported in homology-derived structures of proteins. Beyond the relative entropy for each position in the alignment, ConSSeq also presents the consensus sequence and information about the amino acids, which are predominant at each position of the alignment. ConSSeq is part of the STING Millennium Suite and is implemented as a Java Applet. AVAILABILITY: http://sms.cbi.cnptia.embrapa.br/SMS/STINGm/consseq/, http://trantor.bioc.columbia.edu/SMS/STINGm/consseq/, http://mirrors.rcsb.org//SMS/STINGm/consseq/, http://www.es.embnet.org/SMS/STINGm/consseq/ and http://www.ar.embnet.org/SMS/STINGm/consseq/